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These lectures give an overview of the ongoing application of effective field 
theory (EFT) and renormalization group (RG) concepts and methods to den- 
sity functional theory (DFT) , with special emphasis on the nuclear many-body 
problem. Many of the topics covered are still in their infancy, so rather than 
, a complete review these lectures aim to provide an introduction to the devel- 

^ ' oping literature. 

o ■ 

^ : 1 EFT, RG, DFT for Fermion Many-body Systems 

o ■ 

■ 1.1 Overview of Fermion Many-Body Systems 

P 

There are a wide range of many-body systems featuring fermion degrees 
of freedom. These can be collections of "fundamental" fermions (electrons, 
quarks, ...) or of composites each made of odd number of fermions (e.g., 
protons). Here are some general categories and examples: 



1. Isolated atoms or molecules, which contain electrons interacting via the 
long-range (screened) Coulomb force. 
Jv>( \ 2. Bulk solid-state materials, such as metals, insulators, semiconductors, su- 

5^ ' perconductors, etc. 

5^ I 3. Liquid '^He (a superfluid!). 

4. Cold fermionic atoms in (optical) traps (note that ^Li is a fermion but 
^Li is a boson). 

5. Atomic nuclei. 

6. Neutron stars, which could mean neutron matter or color superconducting 
quark matter. 

DFT has been most widely applied to systems in categories 1 and 2. We will 
focus in these lectures on categories 4 and 5 (and neutrons stars are treated 
in Thomas Schafer's lectures). 
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Fig. 1. Molecule-molecule and nucleon-nucleon potentials compared [2]. 

For our purposes, it won't matter whether the fermions are structureless 
(as far as we know) such as quarks or electrons, or are composites such as 
interacting atoms. Note that an individual atom, studied as a many-body 
system of electrons (with external potential from the nucleus) is a fermion 
many-body system, while a collection of these atoms might be either a boson 
or fermion many-body system [T]. 

If we label the axes appropriately in Fig. [1] (note the many orders of mag- 
nitude difference!), we see qualitative similarities between the central part of 
a conventional nucleon-nucleon (NN) potential and potentials between atoms 
or molecules (a Lcnnard- Jones potential is shown). In particular, there is 
midrange attraction and strong short-range repulsion (or a "hard core"). In 
the atomic case, the attraction is van der Waals in nature from induced po- 
larization (which is why it falls off as 1/r^), and the rapid repulsion is from 
when the electron clouds overlap. In the nuclear case, the long- and mid-range 
attraction is mostly from one- and two-pion exchange. The hard core is often 
described in terms of vector-meson exchange but is generally phenomenolog- 
ical. The potential shown is the central part of the NN interaction; there are 
also important spin dependences and a non-central tensor force [3]. 

What might one expect qualitatively from a many-body system with such 
a potential? Start with the equation of state of an ideal gas PV = nRT 
(n is number of moles). Hard-core means "excluded volume" so V" ^ (V" — 
nb) with b constant. Attraction lowers the pressure on the container, so we 
find: P = ~ i?^- The end result is a van der Waals equation of state: 

(P + ^^){V — nb) — nRT, which has a liquid-gas phase transition. This is 
consistent with nuclei! The hard core keeps particles apart, leading to "short- 
range" correlations in the wave function. They make many-body problems 
difficult but might seem to be essential in saturating the nuclear liquid. Both 
liquid helium and nuclei can be thought of as liquid drops, whose radii scale 
with the number of particles A to the 1/3 power: R ~ vqA^^^. If a hard 
core of radius c is responsible for saturation, one can estimate that 0.55c < 
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Neutron Number ► 

Fig. 2. Overview of low-energy many-body nuclear systems [S]- 

ro < 2.4c [H [S]. Liquid ^He has tq 2.4 A c, but for nuclei ro ~ 1.1 fm ~ 
2.75c. This impHes that nuclear matter is dilute and "delicately" bound, which 
means that an EFT expansion may be particularly useful. Other questions we 
might ask are whether there are common ("universal") features of atomic 
and nuclear systems, and how to we relate the many-body physics to more 
fundamental underlying theories? As we'll see, EFT will also help to address 
these questions. 

In Fig. [21 we show the "big picture" of low-energy nuclear physics, which 
features a specific scientific goal of predicting properties of unstable nuclei 
(the non-black squares), and a general scientific goal, which is connecting the 
whole picture from quantum chromodynamics (QCD) to superheavy nuclei 
in a systematic way with robust predictions. In principle, the nuclear part 
of the problem is simple: given internucleon potentials, just solve the many- 
body Schrodinger equation. This turns out to be feasible only for the smallest 
nuclei. 

Why is this many-body problem so difficult computationally? Let's think 
about the many-body Schrodinger wave function [7]. How can we represent 
the wave function for an A-body nucleus? Consider a *Be (Z = 4 protons, 

= 4 neutrons) wave function with spin, isospin, and space components: 

Yl/) = ^x<tX7-0(R-) , where R are the 3 A spatial coordinates , 

cr,r 

X'y =iih ■■■ Ia (2^ terms) = 256 for A = 8 , 
— nin2 ■ ■ - pA ij^y^ terms) = 70 for *Be . 
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So for ^Be there are 17,920 complex functions in 3A — 3 = 21 dimensions! 
Suppose for a nucleus of size 10 fm you represent this with a mesh spacing 
of 0.5 fm. You would need 10^^ grid points! Obviously we need to drastically 
reduce the necessary degrees of freedom. 

An extreme approximation to the full many-body wave function is the 
Hartree-Fock wave function, which is the best single Slater determinant in a 
variational sense: 

Ii^-hf) =det{0,(x),i = 1---A}, x=(r,a,r). (1) 

The Hartree-Fock energy in the presence of an external potential is [8] 

(i^hfI^I^f) = 




+ E /rfy«cxt(y)|0,(y)r (2) 
1=1 •' 

We determine the by varying with fixed normalization: 

^^((*'HF|i?i'^'HF)-E^^y^y|'^^(y)l') = 0- (3) 

We solve this self-consistently, which is non-trivial because the potential is 
non-local, but drastically simpler than solving for the full wave function. How- 
ever, while Hartree-Fock is a reasonable starting point for atoms, it is not for 
nuclear potentials of the form in Fig. [1] (e.g., note that the Argonne wig "1st 
order" curve in Fig. [5] is not even bound). 

1.2 Density Fiinctional Theory 

An alternative to working with the many-body wave function is density func- 
tional theory (DFT) [Hllinilllli which as the name implies, has fermion den- 
sities as the fundamental "variables" . To date, the dominant application of 
DFT has been to the inhomogeneous electron gas, which means interacting 
point electrons in the static potentials of atomic nuclei. This has led to "ab 
initio" calculations of atoms, molecules, crystals, surfaces, and more 12 . DFT 
is founded on a theorem of Hohenberg and Kohn (HK): There exists an energy 
functional E^^^^ [p] of the density p such that 



(4) 
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Fig. 3. Kohn-Sham DFT for a Vcxt ~ Vno harmonic trap. On the left is the inter- 
acting system and on the right the Kohn-Sham system. The density profile is the 
same in each. 

where ^hk is universal (the same for any external potential Wext); the same 
for i?2 to DNA! This is useful if you can approximate the energy functional. 

The general procedure is to introduce single-particle orbitals and to mini- 
mize the energy functional to obtain the ground-state energy Egs and density 
pgs- This is called Kohn-Sham DFT, and is illustrated schematically in Fig.[3l 
Here the interacting density for A fermions in the external potential ^ho is 
equal (by construction) to the non-interacting density in Vks- Orbitals {V'i(x)} 
in the local potential V^ks([/3],x) are solutions to 

[-Vy2m + VKs{^)]A^e,iP, (5) 

and determine the density 

A 

p(x)==^|V..(x)P (6) 

1=1 

(the sum is over the lowest A states). The magical Kohn-Sham potential 
Vks([p], x) is in turn determined from SEi,^^^ [/3]/(5p(x) (sec below for an exam- 
ple) . Thus the Kohn-Sham orbitals depend on the potential, which depends on 
the density, which depends on the orbitals, so we must solve self-consistently 
(for example, by iterating until convergence). 

DFT for solid-state or molecular systems starts with the HK free energy 
for an inhomogeneous electron gas pLP] : 

i^HK[p(x)] = Tks[p(x)] + ^^l<fx rf3^/ PWP(^0 ^ EM^)] . (7) 

2 J |x-x'| 

Then — ^^cxt — e4> + Wxc with Uxc(x) = SE^d 5p{y.). To calculate the 
Kohn-Sham kinetic energy Tks[/5(x)], find the normalized {tpi^ti} from 

(-^V^ + yKs(x))^.(x) = e,V.(x) (8) 

and, with p(x) = Y.t=i lV'z(x)P, 
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Tks[p(x)] = 5](V',| - |^V?|^.) =J2'^~ [ d'xp{^)VKs{^) . (9) 

In practice, the DFT is usually based on the local density approximation 
(LDA): £'xc[p(x)] w J(Px£^c{p{^)) with £xc{p) fit to a Monte Carlo calcula- 
tion of the uniform electron gas. For example, one parametric formula for the 
energy density is ^lOj 

£M/P - -0.458/r, - 0.0666G(r,/11.4) , (10) 

with 

G(x) = i|(l + x)Mog(l + .T-i)-x2 + ix-i| . (11) 
This is just like a simple Hartree approach with the additional potential: 

d[£^c{p)] 



Wxc(x) 



dp 



(12) 

P=P(x) 



The LDA is improved with the Generalized Gradient Approximation (GGA) , 
such as the van Leeuwen-Baerends GGA [TO] . 

v^4r) = -Pp^/\r) ^^^(J) (^3^ 

with X = \Vp\/ p'^/'^. For these Coulomb systems, Hartree-Fock is generally a 
good starting point, DFT/LDA is better, and DFT/GGA is better stiU. 

There are some concerns, however, about DFT. Here are some quotes from 
the DFT literature that help motivate the application of EFT to DFT: 

From A Chemist's Guide to DFT [T3]: "To many, the success of DFT 
appeared somewhat miraculous, and maybe even unjust and unjusti- 
fied. Unjust in view of the easy achievement of accuracy that was so 
hard to come by in the wave function based methods. And unjustified 
it appeared to those who doubted the soundness of the theoretical 
foundations. " 

From Density Functional Theory 'lOj: "It is important to stress that 
all practical applications of DFT rest on essentially uncontrolled ap- 
proximations, such as the LDA ..." 

From M eta- Generalized Gradient Approximation 14] "Some say that 
'there is no systematic way to construct density functional approxima- 
tions.' But there are more or less systematic ways, and the approach 
taken . . . here is one of the former." 
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Fig. 4. Table of the nuclides. 



Thus, a microscopic, controlled, and systematic approach to DFT would be 
welcome. 

We end this section with a preview of DFT as an effective action [15] . Recall 
ordinary thermodynamics with N particles at T = 0. The thermodynamic 
potential is related to the partition function, with the chemical potential ^ 
acting as a source to change N — (N) , 

f2{fi) ^ -kT In Z{fi) and TV = - ( ^ ) . (14) 

If we invert to find fi{N) and apply a Legendre transform, we obtain 

F{N) ^ n{fi{N)) + fi{N)N . (15) 

This is our (free) energy function of the particle number, which is analogous 
to the DFT energy functional of the density. Indeed, if we generalize to a 
spatially dependent chemical potential J(x), then 

Z(^) — > Z[J(x)] and nN ^ fij i/;^ — > J J(x)V/V'(x) . (16) 

Now Legendre transform from lnZ[J(x)] to /^[p(x)], where p — and 
we have DFT with F simply proportional to the energy functional! 

1.3 DFT for Nuclei: EFT and RG Approaches 

Figure [4] shows the table of the nuclides, labeled by the total numbers of 
protons and neutrons. For example, ^"^Pb is found at the intersection of the 
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Fig. 5. Nuclear matter in perturbation theory for a conventional NN potential 
(Argonne uig) and a low-momentum potential [16j . 



horizontal line labeled "82" (protons) and the vertical line labeled "126" (neu- 
trons). Stable nuclei are in black, known in pink. Here are some basic nuclear 
physics questions (which everyone should know how to answer): 

• Why is the slope of the black region less than a 45 degree angle once it is 
past Z = = 20 or so? 

• How do the binding energies of the nuclei in black compare? (E.g., do they 
vary over a wide range? Do they have a regular pattern?) 

• What happens to the binding energy as you move perpendicular to the 
black line? 

• What is the difference between being unstable and unbound? What are 
the drip lines? 

(See [3j[8] for nuclear physics background.) The nucleosynthesis r-process lives 
largely in "Unexplored Territory." Radioactive beam facilities (existing and 
proposed) address this physics, as well as exotic nuclei such as halo nuclei and 
many other phenomena. As one gets further from stability, the importance 
of pairing grows, which highlights a difference between nuclear many-body 
physics and the physics of Coulomb systems, like atoms and molecules. 

Let's try solving nuclear matter in low-order perturbation theory. One 
expects a minimum in the energy/particle [EjA) versus density (here the 

Fermi momentum fcp oc p /3) at about fcp ~ 1.35 fm"^ and EjA 16MeV. 

The standard Argonne potential [T7] is used in Fig. O with "Brueckner 
ladder" contributions shown order-by-order. First order is Hartree-Fock and 
it is not even bound! The repulsive core means that the series diverges badly. 
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But there is an energy density functional approach that seems to be based 
on Hartree-Fock (HF) , which works quite well throughout Fig. HI This is the 
phenomenological Skyrme HF approach [SI [TSl [111 [50] ■ Recall from our earlier 
Hartree-Fock discussion that solving self-consistently is somewhat tricky be- 
cause the potential is non-local. It would be much simpler if ti(x, y) c>c 5(x— y). 
This is the case with the Skyrme interaction y^Skyrmc _^ ^Skyrmc^ where 

^Skyrme I = to + ^tlCk^ + k'') + tak • k' + iWo{<Ji -|- CJa) • k X k'. This 

motivates the Skyrme energy density functional (for N = Z) [5]: 

+ ^(9ti - 5i2)(Vp)2 - ^WopV ■ J + ^ih - t2)J' , (17) 

where p(x) = J2t and t(x) = X),; (see [^ for the J(x) for- 

mula). We minimize E = Jdx.£[p,T,3] by varying the (normalized) ipi's, 

("^2M^^ + C^(x) + Iwo^p 4^ X ^) M^) = '^i W ' (18) 

with U = |toP + (^^1 + jet2)T -I- • • • and 217%^ 517 + (^^1 + ^*2)p- One 
iterates until the ^^'s and e^'s are self-consistent. 

While phenomenologically successful, there are many questions and pos- 
sible criticisms of Skyrme HF. Typical [e.g., Skill] model parameters are: 
to = -1129, ti = 395, t2 = -95, ts = 14000, Wq = 120 (in units of MeV- 
fm"). These seem large; is there an expansion parameter? Where does 
come from? Why not p^^^l Is this just parameter fitting? A famous quote 
from von Neumann (via Fermi via Dyson) says: "With four parameters I can 
fit an elephant, and with five I can make him wiggle his trunk. " One might 
also say that Skyrme HF is only a mean-field model, which is too simple to 
accommodate NN correlations. How do we improve the approach? Is pairing 
treated correctly? How does Skyrme HF relate to NN (and NNN) forces? And 
so on. There is also the observation that Skyrme functionals works well where 
there is already data, but elsewhere fails to give consistent predictions (and 
the theoretical error bar is unknown). 

Rather than focus on the Skyrme interaction, we consider the Skyrme en- 
ergy functional as an approximate DFT functional. (Note: this is the viewpoint 
of modern practitioners.) Our master plan is to use EFT and renormalization 
group (RG) methods to elevate something close to Skyrme HF to a full DFT 
treatment. We want to use EFT and RG to provide a systematic input poten- 
tial, including three-body potentials, and to generate systematically improved 
energy functionals. At the same time, we want to be able to estimate theoret- 
ical errors, so that extrapolation is under control. 

This is a relatively new and different approach. In Table [1] we summarize 
aspects of the "traditional" approach to the (nuclear) many-body problem 
that are being challenged by the EFT approach. There are many continuing 



10 R.J. Furnstahl 



Table 1. (Nuclear) Many-Body Physics: "Old" Approach 



One Hamiltonian for all problems 
and energy /length scales (not 
QCD!) 


For nuclear structure, protons and 
neutrons with a local potential [21] 
fit to NN data 


Find the "best" potential 


MM i-n-i+ian^-i Q 1 iiri + Vi /r\r\f 1 nri 
IMN pOLcllLld,! WiLil ^ / ClUi ^ 1 Up 

to ~ 300 MeV energy ^ 


Two-body data may be sufficient; 
many-body forces as last resort 


Let phenomenology dictate 
whether three-body forces are 
needed (answer: yes! [23)) 


Avoid (hide) divergences 


Add "form factor" to suppress 
high-energy intermediate states; 
don't consider cutoff 
dependence [24] 


Choose approximations by "art" 


Use physical arguments (often 
hand waving) to justify the subset 
of diagrams used [25] 



successes of conventional many-body approaches. The idea is not to prove 
standard methods wrong but to highlight where new insight can be provided. 
For each "old" item in this table (see endnotes for further explanations), we'll 
have a "new" perspective from EFT (see Table [2]). 

1.4 EFT Analogies 

From an effective field theory we desire systematic calculations with error 
estimates. We want them to be robust and model independent, which will 
enable reliable extrapolation. To help understand how this is accomplished, 
we can explore useful analogies between EFT and sophisticated numerical 
analysis. 

• Naive error analysis: pick a method and reduce the mesh size (e.g., increase 
grid points) until the error is "acceptable" . Sophisticated error analysis: 
understand scaling and sources of error (e.g., algorithm vs. round-off er- 
rors). Does it work as well as it should? 

• Representation dependence means that not all ways of calculating are 
equally effective! 

• Reliable extrapolation requires completeness of an expansion basis. An 
EFT lagrangian provides the analog of a complete basis. 

Note: quantum mechanics makes EFT trickier than "classical" numerical anal- 
ysis (see Lepage's lectures)! 

Consider the numerical calculation of an integral using equal-spaced inte- 
gration rules of increasing sophistication. How do the numerical errors behave? 



EFT for DFT 11 



Numerical Integration 



2 

+ trapezoid rule 0(h ) 

4 

^ Simpson's rule 0(h ) 

6, 



10' 

2 10'"^ * Milne's rule 0(h'') 
o + 

> -12 + 

•■= 10 - ^ + 

10-'^- 



10-' 



10-' 



10"^ 



10 10 
mesh size h 



* * 

I 

■3 



X * 



10" 



10-' 



Fig. 6. Error plots in numerical integration. 
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;. 7. The representation can make a difference: errors for (|19|l vs. (|20|) 



Log-log plots of the relative error against a parameter such as the mesh size 
are very helpful; a straight line indicates a power law and the exponent is read 
off from the slope. An example is shown in Fig. [51 Similar plots can be made 
for order-by-order EFT calculations, as described in Peter Lepage's lectures. 
Just like computer math is not equivalent to ordinary math, EFT is not the 
same as a theory applicable at all energies. It breaks down at high resolution, 
as the numerical calculations break down and degrade because of round-off 
errors at small mesh sizes. (Note: Don't carry this analogy to extremes!) Next 
consider an elliptic integral: 

/ y/{l~x^){2-x) dx (19) 
Jo 

with errors plotted in Fig. [71 Something is wrong; the errors do not behave as 
expected. However, after a simple transformation: 
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/ sin^ y y/2 — cos y dy . (20) 
Jo 

As seen in Fig. [71 the transformation makes a big difference! ff you have 
freedom to change the representation, you can make a big difference in the 
abihty to calculate accurately and easily. This is a moral we will apply when 
using EFT and RG methods to many-body problems. 

1.5 Principles of Effective Low-Energy Theories 

There are some basic physics principles underlying any low-energy effective 
model or theory. A high-energy, short-wavelength probe sees detail down to 
scales comparable to the wavelength. Thus, high-energy electron scattering at 
Jefferson Lab resolves the quark substructure of protons and neutrons in a 
nucleus. But at lower energies, details are not resolved, and one can replace 
short-distance structure with something simpler, as in a multipole expansion 
of a complicated charge or current distribution. So it is not necessary to do full 
QCD to do strong interaction physics at low energies, we can replace quarks 
and gluons by neutrons and protons (and pions). Chiral effective field theory 
is a systematic approach to carrying out this program using a local lagrangian 
framework. 

It is not obvious that working at low resolution will work in quantum me- 
chanics as it does for pixels or point dots or the classical multipole expansion, 
because virtual states can have high energies that are not, in reality, simple. 
But renormalization theory says it can be done! (See Lepage's lectures and 
[26l [27].) Note that this doesn't mean that we are insensitive to all short- 
distance details, only that their effects at low energies can be accounted for 
in a simple way. 

We can use the possibility of changing the resolution scale to change the 
"perturbativeness" of nuclei. There are several sources of nonperturbative 
physics for nucleon-nucleon interactions: 

1. Strong short-range repulsion ("hard core"). 

2. Iterated tensor {S12) interactions (e.g., from pion exchange). 

3. Near zero-energy bound states. 

In Coulomb DFT, Hartree-Fock gives the dominate contribution to the en- 
ergy, and correlations are small corrections. This may be why DFT works. In 
contrast, for NN interactions, correlations ^ HF; does this mean DFT fails?? 
However, the first two sources depend on the resolution (i.e., the cutoff of 
high-energy physics), and the third one is affected by Pauli blocking. Thus we 
might use the freedom of low-energy theories to simplify calculations. 

We can see the impact of different resolutions on the deuteron wave func- 
tion in Fig. [H The repulsive core leads to short-distance suppression and im- 
portant high-momentum (small A) components. This makes the many-body 
problem complicated. In contrast, potentials evolved by renormalization group 
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Fig. 8. The deuteron probability density at different resolutions (as indicated by 
the sharp momentum cutoff A). 
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(RG) methods to low momentum, generically called Mowfc [30l [311 [32] , have 
much simpler wave functions! (See Andreas Nogga's lectures for more de- 
tail on such potentials.) We note that chiral EFT potentials are naturally 
low-momentum potentials, but lowering their cutoffs further is generally ad- 
vantageous. The consequence of lower resolution for variational calculations 
is illustrated in Fig. [9| [28l [29] ■ Note that the improvement for the deuteron 
comes with smooth (exponential) cutoffs '33|, which is another example of 
how the representation can make a difference. 

These observations carry over to nuclear matter as well, as seen in Fig. [10] 
(although the effect with lowered cutoff in the important '^Si channel is less 
dramatic). In medium, the phase space in the pp-channel is strongly sup- 
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Fig. 10. Two-body correlations in nuclear matter. 
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(In addition, the potential itself gets smaller in magnitude in the integration 
region.) This tames the hard core, tensor force, and the bound state. If we 
return to Fig. \S\ we see the consequence, which is very rapid convergence 
(by 2nd order) of the in- medium T- matrix for Mowfe [Hj. But there is no 
saturation in sight! 

There were active attempts to transform away hard cores and soften the 
tensor interaction in the late sixties and early seventies [SH [3S1 . But the 
requiem for soft potentials was given by Bethe (1971): 

"Very soft potentials must be excluded because they do not give satu- 
ration; they give too much binding and too high density. In particular, 
a substantial tensor force is required." 

The next 30-1- years were spent trying to solve accurately with "hard" po- 
tentials. But the story is not complete: three- nucleon forces (3NF)! When 
they are added consistently, we have the advantages of soft potentials while 
answering Bethe's criticism. 

Ideally we would start with chiral NN -I- 3NF EFT interactions and then 
evolve downward in A. What is possible now is to run the NN and fit 3NF 
EFT at each A [37j . The consequence is shown in Fig. [TT] There is saturation 
even at the Hartree-Fock level, which is now driven by the three-body force. 
At second order, the cutoff dependence is greatly reduced and the minimum 
moves closer to the empirical point. One might worry that the three-body force 
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is unnaturally large, but it is consistent with EFT power counting. Excellent 
results are also found for neutron matter. While there is much to be done, 
these results are very encouraging and motivate a microscopic DFT for nuclei. 

1.6 Summary 

In summary, there is a new attitude for many-body physics inspired by effec- 
tive field theory. Bethe wrote about the nuclear case: 

"The theory must be such that it can deal with any nucleon-nucleon 
(NN) force, including hard or 'soft' core, tensor forces, and other com- 
plications. It ought not to be necessary to tailor the NN force for the 
sake of making the computation of nuclear matter (or finite nuclei) 
easier, but the force should be chosen on the basis of NN experiments 
(and possibly subsidiary experimental evidence, like the binding en- 
ergy of H^)." 

The new attitude is to instead seek to make the problem easier. It's like the 
old vaudeville joke about a doctor and his patient: 
Patient: Doctor, doctor, it hurts when I do this! 
Doctor: Then don't do that. 

We also follow Weinberg's Third Law of Progress in Theoretical Physics: "You 
may use any degrees of freedom you like to describe a physical system, but 
if you use the wrong ones, you'll be sorry!" We conclude with a new table 
(Table [2]) of many-body physics, contrasting the old and the new approaches. 
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Table 2. (Nuclear) Many-Body Physics: "Old" vs. "New" 



One Hamiltonian for all problems 
and energy /length scales 


Infinite # of low-energy potentials; 
different resolutions =>different 
dof's and Hamiltonians 


Find the "best" potential 


There is no best potential 
=>use a convenient one! 


Two-body data may be sufficient; 


Many-body data needed and 


many-body forces as last resort 


many-body forces inevitable 


Avoid divergences 


Exploit divergences (cutoff 
dependence as tool) 


Choose diagrams by "art" 


Power counting determines 
diagrams and truncation error 



2 EFT/DFT for Dilute Fermi Systems 
2.1 Thermodynamics Interpretation of DFT 

As an analogy, consider a system of spins Si on a lattice with interaction 
strength g j38| . The partition function has all the information about the energy 
and magnetization of the system: 



The magnetization M is 



=lTr 
Z 



(23) 

(24) 
(25) 



Now add an external magnetic probe source H. The source adjusts the spin 
configurations near the ground state, 



Variations of the source yield the magnetization 



H 



dH 



(26) 



(27) 



where F[H] is the Helmholtz free energy. For the ground state, we set H = 
(or equal to a real external source) at the end. 
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Now if we find H[AI] by inverting AI[H], 



H 



dF[H] 
dH 



(28) 



we can Legendre transform to the Gibbs free energy 

r[M] ^ F[H] + H M . (29) 
Then the ground-state magnetization il/gg follows by minimizing _r[Af]: 

dr[M] dr[M] 



H = 



dM 



DM 



= 



(30) 



Thus, we have a function of M (or functional if H is inhomogeneous) that is 
minimized at the ground-state free energy and magnetization. 

DFT has an analogous structure as an effective action [TS]. An effective 
action is generically the Legendre transform of a generating functional with an 
external source (or sources). For DFT, we use a source to adjust the density 
instead of the magnetization. The partition function in the presence of J{x) 
coupled to the density is (we'll use a schematic notation here): 



Z[J]=e-^'W ~Tre-^(^+'^^) 



V 



[i>^]V[i>]e-I^^+''^'^K (31) 



The density p{x) in the presence of J{x) is (keep in mind that we want J ~ 
eventually). 

After inverting to find J[p], we Legendre transform from J to p: 

6r[p] 



r[p] = W[J] - j J P and J{x) 



Sp{x) 



(33) 



Now consider the partition function in the zero-temperature limit of a 
Hamiltonian with time- independent source J(x) [551 140j : 



(34) 



H{J) = H + J Jip^ip ■ 
If the ground state is isolated (and bounded from below) 

)(0|+O(e-^(^i-^°)) 



As /? ^ cx), Z[J] yields the ground state of H{J) with energy Eq{J): 
Eo{J)= \im -^\ogZ[J] = ^W[J]. 



(35) 



(36) 
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Substitute and separate out the pieces: 

Eo{J) = {H{J))j = {H)j + [j{^H)j = {H).j + / Jp(J) . (37) 



Rearranging, the expectation value of in the ground state generated by J[p] 
is 

{H)j = E^{J)~ [ Jp^\r[p]. (38) 



/3 



Let's put it all together: 



\r[p\ = {H)j'-^E, and J(.) ^np] J^o 5r[p] 



P ' 5p{x) Sp{x) 



= . (39) 



So for static p(x), r[p] is proportional to the DFT energy functional Fhk' 
Furthermore, the true ground state (with J = 0) is a variational minimum, so 
additional sources should be better than just one source coupling to the den- 
sity (as we'll consider below). The simple, universal dependence on external 
potential follows directly in this formalism: 

r^p] = Wy[J]-Jjp^W,^o[J+v]-J[{J+v)~v]p^r,=o[p]+ Jvp . (40) 

[Note: the functionals will change with resolution or field redefinitions; only 
stationary points are observables.] 

There are a number of paths to the DFT effective action: 

1. Follow the usual Coulomb Kohn-Sham DFT by calculating the uniform 
system as function of density, which yields an LDA ("local density ap- 
proximation" ) functional and a standard Kohn-Sham procedure. Improve 
the functional with a semi-empirical gradient expansion. 

2. Derive the functional with an RG approach [41] . 

3. Use the auxiliary field method [421 [43]. Couple V'^V' to an auxiliary field 
ip, and eliminate all or part of (V'^V')^- Add a source term jLp and perform 
a loop expansion about the expectation value (fi = {(p). A Kohn-Sham ver- 
sion uses the freedom of the expansion to require the density be unchanged 
at each order. 

4. The inversion method [lH [45l [46l [47] yields a systematic Kohn-Sham 
DFT, based on an order-by-order expansion. For example, we can apply 
the EFT power counting for a dilute system. 

We'll expand here on the last path. 



2.2 EFT for Dilute Fermi Systems 

We consider first one of the simplest many-body systems: a collection of "hard 
spheres," which means that the potential is infinitely repulsive at a separation 
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sin(kr+5) 






Fig. 12. Hard-sphere phase shifts and scales at finite density. 



R of the fermions. Since the potential is zero outside of R and the wave func- 
tion must vanish in the interior of the potential (so that the energy is finite), 
we can trivially write down the 5- wave scattering solution for momentum k: it 
is just a sine function shifted by kR from the origin (see Fig.ll2p. Our problem 
will be to find the energy per particle (and other observables) of a system of 
particles interacting with this potential at T = 0, given the density. 

Let's do a quick review of scattering. (More details on scattering at this 
level can be found in practically any first-year graduate quantum mechanics 
text. For a more specialized but very readable account of nonrelativistic scat- 
tering, check out "Scattering Theory" by Taylor.) Consider relative motion 
with total momentum P — 0: 



„ikr 

fik,e)— 

r 



(41) 



where k'^ ^ k'"^ — MEk and cos6' = k ■ k' . The differential cross section is 
da/dfl — |/(fc,6')p. For a central potential, we use partial waves: 



where l48l 



/(fc,0)=^(2/ + l)/Kfc)fl(cose), 
I 

f,tSi(k) sin 6i{k) _ 1 



k k cot 6i (fc) — ik 

and the 5- wave phase shift is defined by 

'"o('') sin[fcr -|- SQ{k)] =^ ^oik) = —kR. for hard sphere . 



(42) 



(43) 



(44) 



Note: we can do a partial wave expansion even if the potential is not central 
(as in the nuclear case!); it merely means that different Vs will mix. The more 
important question is how many total Vs do we need to include to ensure 
convergence. 

As first shown by Schwinger, fc'+^cot5;(fc) has a power series expansion. 
For I = 0, 
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fc cot (5o = - — + irofc2 ^Prlk'^ + ... , (45) 
ao 2 

which defines the scattering length ao and the effective range tq. While tq ~ i?, 
the range of the potential, ao can be anything; if ao ^ R, it is called "natural" . 
The other case |ao| 3> R (unnatural) is particularly interesting; it is the case 
for nucleon-nucleon interactions and can be studied in detail with cold atoms. 
The effective range expansion for hard sphere scattering is: 

kcot{-kR) = -^ + \Rk'^ + ■■■ =^ aQ = R ro = 2i?/3 , (46) 
R 3 

so the low-energy effective theory is natural, is Schwinger first derived the 
effective range expansion back in the 1940's and then Bethe showed an easy 
way to derive (and understand) it. The implicit assumption here is that the 
potential is short-ranged; that is, it falls off sufficiently rapidly with distance. 
This is certainly satisfied by any potential that actually vanishes beyond a 
certain distance. Long-range potentials like the Coulomb potential must be 
treated differently (but a Yukawa potential with finite range is ok) . 

So now consider the EFT for a natural, short-ranged interaction [l^. A 
simple, general interaction is a sum of delta functions and derivatives of delta 
functions. In momentum space, 

(kl^eftlk') =Co + ic2(k2 + k'') + k • k' -f • • • (47) 

Or, we construct the effective lagrangian Cdt from the most general local 
(contact) interactions: 



2M 



^-^(V't7^)' + §[(V'V')^('^v'^)+ h.c.; 



+ ^(V^vV')t • (VvV-) - ^(^V)' + . . . (48) 

8 6 

Dimensional analysis (with a bit of additional insight to give us the 47r's) 
implies 

C,.-^R'^+\ D,,^^R'^+\ (49) 

which will enable us to make quantitative power-counting estimates. 

The ingredients for an effective field theory are nicely summarized in the 
"Crossing the Border" review [5D]: 

f. Use the most general C with low-energy degrees- of- freedom consistent with 
global and local symmetries of underlying theory. Here, 
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P/2 + k P/2 + k' 




P/2-k P/2-k' 

-8(k'|FEFT|k) 





+ 



k2 _^ jj./^ 

—iCo —iC2 — ■iC2k-k' 



-iDr. 



Fig. 13. Feynman rules in free space [49] . 

2. Declare a regularization and renormalization scheme. For a natural oq, 
using dimensional regularization and minimal subtraction is particularly 
convenient and efficient. 

3. Establish a well-defined power counting, which means identifying small 
expansion parameters, typically using the separation of scales. Here, kp/ A 
with A ~ which implies fepag, kpro, etc. are expansion parameters. 
In the end, this will be manifest in the energy density: 



c F 



I + ^{kpao) + - 2 In 2)(fcFao)' 



(51) 



The Feynman rules for the EFT lagrangian are summarized in Fig. flSl [49] . 
We need to reproduce /o(fc) in perturbation theory (the Born series): 



/o(fc) (X ao - ia^k - (ag - a^ro/2)k^ + ©(fc-^a^J) 
The leading potential V^prp(x) = Co(5(x) or 



(52) 




Co 



(53) 



Thus, choosing Cq (x oq gets the first term. Next is (k|V^GoV^|k'): 

d^q 1 




CnM 



{2n f k^-q^+ i, 



■Cn 



(54) 



This is a linear divergence. If the integral is cutoff at Ac, we can absorb the 
linear dependence on Ac into Cq, but we'll have all powers of fc^: 



d\ 



1 



(27r)3 k^-q^ + it 



27r2 4^ 



(55) 
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A more efficient scheme is dimensional regularization with minimal subtrac- 
tion (DR/MS), which implies only one power of k survives: 



/ 



d^q 1 D^3 ik 

(56) 



The diagrammatic power counting with DR/MS is very simple, with each 
loop adding a power of k: 



P/2 + k P/2 + k' 






P/2-k P/2-k' 

iT(k,cos9) -iCa -—{Cafk 

47r 






0(k^ 



iCnk^ -iCi,k^cos9 



After matching to the scattering amplitude, 

recovers the effective range expansion ordcr-by-ordcr with diagrams; 

^ (ao - ialp - alp^ + agrop^ H ) , (58) 

with one power of k per diagram and natural coefficients, so we can estimate 
truncation errors from simple dimensional analysis. 

2.3 Apply at Finite Density 

Consider a noninteracting Fermi sea of particles at T = 0. Put the system 
in a large box {y = L^) with periodic boundary conditions and spin-isospin 
degeneracy v (e.g., for nuclei, i/ = 4). Fill momentum states up to Fermi 
momentum k-p, so that 

fep fcp 
k k 



We can evaluate the sums using 

= F(kA _ 

L L 



I 



27r 27r 

F{k)dk^Y.F{ki)Aki = Y,F{ki) — Ani = —Y,F{h) . (60) 
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In one dimension (try finding the 1-D analogs of the foUowing results!), 

while in three dimensions: 

N = v- — -- / (Pk= — =^ p= — = — ^; — = E-p . (62) 

(27r)3 J 67r2 ^ V Qtt^ V ^ 2M ^ ^ ' 

The volume/particle VjN = 1/ p ^ l/fc|, so the spacing ~ l/fcp, as implied 
by Fig. m 

We find the energy density by summing over the Fermi sea. In leading 
order, we found V^plp(x) = Co5(x), so that V^plp(k, k') = Cq, and 

OO fLO = ^K^-l)(^£l^ (xaofc^ (63) 
At the next order, we get a linear divergence again: 

} ^--IIfA- '''' 

The same renormalization fixes it! 

1 0-^3 f^" 1 2.7 

12 2 I 12 2 I 71 2 ^ / T2 2 '^0% ■ 

« " 9 Jo « - 9 Jo - 9 Jo «^ - Q 

We also note that particles — > holes through the renormalization. 

The Feynman rules for the energy density £ at T = 0, which is the sum 
of HugenhoUz diagrams |38| (closed, connected Feynman diagrams with sym- 
metry factors) with the same vertices as free space (and the same renormal- 
ization!), are: 

1. Each line is assigned conserved k = (fco, k) and [ujk = k'^/2M], 

■^,7, T ( Oik-kp) d{kF-k) \ 
\ fco - Wfc + ie ko - LJk - te J 

s s 

5[js + SasSfj^) (if spin- independent). 

3. After spin summations, 5aa — > in every closed fermion loop. 

4. Integrate /d'*fc/(27r)'* with e'''«°* for tadpoles 

5. The symmetry factor i/('^n;=2''(^')'^) counts vertex permutations and 
equivalent /-tuples of lines (see [H] for examples). 

These Feynman rules in turn lead to power counting rules: 
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1. for every propagator (line): M/I^p 

2. for every loop integration: A^/M 

3. for every n-body vertex with 2i derivatives: Aj^/Afyl^*+'^"~^ 
Thus, a diagram with VJ^ n-body vertices scales as (kp)^ only: 

oc oc 

/3 = 5 + ^^(3n + 2z-5)y2". 



(67) 



For example, at leading order, 



OO 



/3 = 5 + (3-2 + 2-0-5)-l = 6=^C'(4) , (68) 



and at next-to-leading order. 



yQ2^2=4>/3 = 5+(3-2 + 2- 0-5)-2 = 7=^ 0(4) • (69) 



We emphasize that Pauli blocking doesn't change the free-space ultraviolet 
(short distance) renormalization, since the density is a long-distance effect. As 
noted before, particles become holes: 



(70) 



The power counting is exceptionally clean, with a separation of vertex fac- 
tors cx flQ, ro, . . . and a dimensionless geometric integral times fcp, with each 
diagram contributing to exactly one order in the expansion. This is a system- 
atic expansion: the ratio of successive terms is ~ kpR, so you can estimate 
excluded contributions. 

The full result for the energy density through O(fcp) is [49] : 



E 
V 



'2M 



^ + - l)|^(fcFao) + 1)^(11 - 21n2)(fcFao)' 



+ 1)(0.076 + 0.057(1/ - 3)) (kpaof + (i^ - 1) — (A:Fro)(/cFao)^ 



+ {u + l) — {kFapf + 

OTT 



(71) 



This looks like a power series in fcp, but it's not! There are new logarithmic 
divergences in 3-3 scattering, in these diagrams: 





cx(Co)^ln(fcMe) 



(72) 
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Changes in must be absorbed by the 3-body couphng Dq{Ac), so [5T] 

Dq{Ac) oc {Cof \n{aoAc) + const. (73) 
Then requiring the sum to be independent of Ac, 



d 

dAr. 









(74) 



fixes the coefficient! This imphes for the energy density, 




ca [v ~ 2)(z^ — 1) {k-pao)'^ In(fcFao) 



(75) 

without actually carrying out the calculation. Similar analyses can identify 
the higher logarithmic terms in the expansion of the energy density [51] [49j . 

Divergences indicate sensitivity to short-distance behavior. The cutoff A^ 
here serves as a resolution scale; as we increase A^., we see more of the short- 
distance details. Observables (such as scattering amplitudes) must not change 
when Ac changes, so they must be absorbed in a coupling. But it can't be a 
coupling from 2-2 scattering, because we already took care of all the diver- 
gences there. So there must be a point-like three-body force included, whose 
coupling Do can absorb the changes. 

Let's summarize the dilute Fermi system with natural ao. We find that 
the many-body energy density is perturbative (but not analytic) in fepoo: 
and is efficiently reproduced by the EFT approach. Power counting gives us 
error estimates from omitted diagrams. Three-body forces are inevitable in 
the low-energy effective theory and not unique; they depend on the two-body 
potential. The case of a natural scattering length is under control for a uniform 
system, but what if the scattering length is not natural? We'll come back to 
that situation in the last lecture. First we consider a non-uniform system: a 
finite number of fermions in a trap, which takes us back to DFT. 



2.4 DFT via EFT [521, [53l [54] 

Return to the thermodynamic version of DFT through the effective action 
and ask: What can EFT do for DFT? We can construct the effective action 
as a path integral by finding W{J] order-by-order in an EFT expansion. For a 
dilute short-range system, this means the same diagrams as before, but now 
the propagators (lines) are in the background field J(x): 

G",(x,x»=^^„(x)V':(x') 

a 

where V'a(x) satisfies: 



0{ea - ef) 



(ep - £«) 



irj 



IT] 



(76) 
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V2 

[-2M ^ ^'^'^t W - ^W] '/'^(x) = eQ'(/'a(x) . (77) 

Applying this to the leading-order (LO) contribution Wi[J], which is Hartree- 
Fock, yields 

W,[J] = \.{u-l)C, jd'^ G",(x,x;c.)G",(x,x;c.') 

= -Y^C, j d^^[pj{^)Y , (78) 

where p,j{x.) = jyj |V'a(x)p. Expressions for the other Wi's proceed directly 
from the Feynman rules using the new propagator. 

Given W[J] as an EFT expansion, how do we perform the Legendre trans- 
formation, 

r[p] - w[j] - jjp. (79) 

in a systematic way? The EFT power counting gives us a means to invert to 
find J[p]. In particular, the "inversion method" provides an order-by-order 
inversion from W[J] to r[p\ [44l |45l |46] . It proceeds by decomposing J{x) = 
Jo{x) + Jho{x) + JNLo(a^) + • ■ • with two conditions on Jq: 



■ (80) 

P=Pg3 



We are using the freedom to split J into Jq and the rest, in the same way 
that one adds and subtracts a single-particle potential [/ to a Hamiltonian: 
H = T + V = [T + U) + [V -U) and then uses the freedom to choose U to 
improve many-body convergence. In our case, we choose Jq so that there are 
no corrections to the zeroth order density at each order in the expansion. The 
interpretation is that Jo is the external potential that yields for a noninter- 
acting system the exact density. This is the Kohn-Sham potential! The two 
conditions involving Jq imply a self-consistent procedure. 

The inversion method for effective action DFT [44l |45l [46] is an order- by- 
order matching in a counting parameter A (e.g., an EFT expansion): 

diagrams =^ W[J, A] ^ Wo[J] + \Wi [J] + A^PFa [J] + • • • (81) 

assume =^ J[p, A] = Jo[p] + XJi[p\ + A^/2[p] H (82) 

derive =^ r[p, A] = Tq [p] + AA [p] + X^T^ [p] + • • • (83) 

We start with the exact expressions for F and p [note: (3 and T — 1 here], 

r[p\ = W[J] - j d^x J{x)p{x) , (84) 



(85) 
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Then plug in each of the expansions, with p treated as order unity. Zeroth 
order is the noninteracting system with potential Ja{x)^ 



p{x) 



SJo{x) 



(86) 



which is the Kohn-Sham system with the exact density! (Note: Jq = Vks here.) 
To evaluate Wo[Jo], we introduce the orbitals from ([77]) . which diagonalize Wq, 
so that it yields a sum of e^'s for the occupied states. Then we find Jq for the 
ground state via a self-consistency loop: 



Jo ^ T^i ^ A ^ Ji W2 



Mx)^J2 



i>0 



sr,[p] 

6p{x) 



(87) 



which is the second of our two conditions on Jq. 
We note that the Kohn-Sham potential is local: 



Jo(x) - 



<5p(x) 



in stark contrast to the non-local and state-dependent self-energy 17* (x, x'; w). 
Evaluating the functional derivatives is easiest if F is approximated so that 
the dependence on the density is explicit, as with the LDA or DME (see 
below). Otherwise we need to use a chain rule with the "inverse density- 
density correlator" [45] 



Jo(R) 



(^Tiiit[p] 
<5p(R) 



■I U,7„(y)j 5J„(y) 




There are new Feynman rules for /]nt for evaluating such diagrams [45]. (A 
related approach is the OEP method [TTIISSIISSIIST] .) 

In constructing the diagrams for W[J] and new diagrams for r[p\ order by 
order in the expansion (e.g., EFT power counting), the source Jo(x) is now 
the background field (rather than the full J(x)). Propagators (lines) in the 
background field Jo(x) are 



G^s(x,x';c.)=5 

c 

where V'a(x) satisfies: 



I Va(x)V'*(x') 



ep) , 0{fF - ea) 



+ 



irj 



V2 

+ i'(x) - Jo(x)] ?/'a(x) = e„?Aa(x) 



(89) 



(90) 



For example, if we apply this prescription to the short-range LO contribution 
(i.e., Hartree-Fock), we obtain 
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1(^-1) 



Co hf^[pj,X^)Y 



where 



(91) 



(92) 



Let us construct the T = local density approximation (LDA). In a uni- 
form system, each line is a non-interacting propagator. The energy density in 
the uniform system evaluates to: 



E 



V ' 2M 



- + 1) — (fcpao) + iiy~ - 21n2)(fcFao) 



357r2 



+ (i/ - 1) (0.076 + 0.057(1/ - 3)) (kpaof + {i^ - 1) — (fcFro)(/cFao)^ 



+ {„ + l) — {kFap)^ + 



(93) 



with fcp = {6tt'^ p/v)-^^^ . Using this relation to replace fcp everywhere by p{x), 
we directly obtain the LDA expression for r[p], 



r[p] 



Tks (x)-^! ^ [p(x)]^ + d,^ [p(x)]^/3 



+ d2ag[p(x)]«/3 + rf3a§ro[p(x)]8/3 
+ d4 4[p(x)]«/3 + . 



(94) 



The Kohn-Sham Jq according to the EFT expansion follows immediately in 
the LDA from 



Jo(x) 



^^^W-cr||[p(x)]V^-c.a3[,(.)]5/3 



-C3agro[p(x)]5/3-c4a3[p(x)f/3 + ... . (95) 

(Finding the {(ii}'s and {c^j's is left as an exercise for the reader.) 
Given ([M)) and ([55]) . the iteration procedure is: 

1. Guess an initial density profile p{r) (e.g., the Thomas- Fermi density). 

2. Evaluate the local single-particle potential VKs(r) = Uext(^) ~ '^o('')- 
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4 



= (exact) 




Kohn-ShamLO 

Kohn-Sham NLO (LDA) 

Kohn-Sham NNLO (LDA) 



E/A <kpa^> <r > 
6.75 -0.52 2.60 



1/2 



5.98 -0.58 2.35 



-0.55 2.47 - 



6.23 -0.55 2.46 



0, 







2 



3 



4 



5 



r/b 



Fig. 14. Density profile of a dilute system of fermions in a trap [52] , 

3. Find the wave functions and energies {ipa, Cq} of the lowest A states (in- 
cluding degeneracies) by solving: 



4. Compute a new density pir) — J2a=i IV''a(x)p. Other observables are 
simple functionals of {V'q, 

5. Repeat 2.-4. until changes are small ("self-consistent") 

This sounds like a simple Hartree calculation! Results at different EFT orders 
for a dilute Fermi gas in a harmonic oscillator trap is given in Fig. [M) Note 
the systematic progression from order to order. 

An important consequence of the systematic EFT approach is that we can 
also estimate individual terms in energy functionals. If we scale contributions 
to the energy per particle according to the average density or (fcp), we can 
make estimates [52l[53]. This is shown in Fig. [15] for both the dilute trapped 
fermions, which is under complete control, and for phenomenological energy 
functionals for nuclei, to which a postulated QCD power counting is applied. 
In both cases, the estimates agree well with the actual numbers (sometimes 
overestimating the contribution because of accidental cancellations), which 
means that truncation errors are understood. 

Conventional DFT is one example of using effective actions, which feature 
sources coupled to composite operators. It's possible that for some apphcations 
a different type of effective action may be better. There are many outstanding 
questions from the present discussion, particularly as we try to adapt it to 
real nuclei. We'll address some of them in the next lecture. 
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Fig. 15. Estimates for energy functionals for a dilute fermions in a harmonic trap 
(left) and for phenomenological energy functionals for rmclei (right). 



3 Refinements: Toward EFT/DFT for Nuclei 

Let's enumerate some questions about DFT and nuclear structure. 

• How is Kohn-Sham DFT more than mean field? That is, where are the 
approximations and how do we truncate? How do we include long-range 
effects (correlations)? 

• What can you calculate in a DFT approach? Can we calculate single- 
particle properties? Or excited states? 

• How does pairing work in DFT? Can we (should we) decouple pp and phi 
Are higher-order contributions important? 

• The Skyrme functional depends on multiple densities: p(x), t(x), and J(x); 
how does that work? 

• What about broken symmetries that arise with self-bound systems? (trans- 
lation, rotation, . . . ) 

• How do we connect to the free, microscopic NN- • -N interactions? Can we 
use chiral EFT or low- momentum intcractions/RG? 

We'll explore some answers to these questions (and note which ones are open) 
in this lecture. 

Consider Kohn-Sham DFT compared to the Thomas-Fermi energy func- 
tional, for which the entire functional is treated in the local density approxi- 
mation (LDA). In Kohn-Sham DFT, treating kinetic energy non-locally leads 
to the shell structure of electrons in atoms and in trapped atoms, as seen in 
Fig. [m This motivates going even further beyond the LDA. 

As a simple step beyond Kohn-Sham LDA, we consider functionals of 
the kinetic energy density in addition to the usual fermion density. The phe- 
nomenological Skyrme E is a. functional of p and r = (Vt/i^ • Vi/") (and J): 

E[p, rj]^l d^x I + ^top' + ^*3P'+" + ^(3ti + 5<2)PT 
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Kohn-Shara NNLO (LDA) 

- - - Thomas-Fermi NNLO (LDA) 




E/A <kpa^> 
2.849 -0.358 
2.785 -0.362 



Fig. 16. Thomas- Fermi vs. DFT for atoms [10) (left) and trapped fermions (right). 



+ ^(9ti - 5i2)(Vp)2 - -^W^pV ■ J + ^(^1 - t2)J'| . (97) 

To do this in DFT/EFT, add to the Lagrangian r]{yi) V^^ and generalize 
our Legendre transformation and inversion to r[p, r], 



r[p, t] - W[J, 77] - / J{x)p{x) - / ?7(x)t(x) 



Now there are two Kohn-Sham potentials 

sr-.nt [p, T 



^o(x) 



and r/o (x) 



^^int [P,t] 



(98) 



(99) 



(5p(x) (5r(x) 
The Kohn-Sham equation defines l/A/*(x) = 1/M — 2ryo(x): 

(-V • ^^j^^V + i;ext(x) - Jo(x)) (/)a(x) = e„ (/.„(x) . (100) 

A simple first application is to evaluate Hartree-Fock diagrams including 
the quadratic gradient terms [S3j. Consider the HF "bowtie diagrams" 



CX3 + CK3 



that have vertices with derivatives: 
, C2 



^cft 



16 



The energy density in Kohn-Sham LDA is 



^int [P] 



C2 r3 /67r 



2 \ 2/3 



8/3 



+ 



3 /67r_^ 
5 W 



2\ 2/3 



r,8/3 



(101) 



(102) 



while the energy density in Kohn-Sham with t {v = 2) is 
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LO NLO LDA pt 10*Vp LO NLO LDA pt 10*Vp 

X -NNLO T -NNLO 

Fig. 17. Estimates of terms in the energy functional, including those with gradients, 
compared to actual values [53) . 




fint [P, r] = ... + ^[pr+ \{Vpf] + ^ [pr - \{V pf] +... (103) 

We find that power counting estimates for terms in the energy functional also 
work with gradient terms (see Fig. [TT]). 

Now let's compare the dilute fermion functional to the phenomenological 
Skyrme functional. The Skyrme energy density functional (for N — Z) is given 
in (j97p while the corresponding dilute energy density functional for v ^ A (and 

Kixtcrnal — 0) is 

eIp.t.j] = Jd'iS^ + ^Cw'' + -^{Xh + 5Ci)i,T+-^{9C;-5c:MVpf 
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Fig. 19. Comparing effective masses and single-particle spectra from energy func- 
tional of p only and p, r [53) . 



They have the same terms after the association ti ^ Ci, except that the 
Skyrme functional is missing non-analytic terms, the three-body contribution, 
and other features. We note after matching to empirical Skyrme coefficients 
that the "effective" scattering length from Co is oq « — 2-3fm, but that 
jfepapl, I^F^'ol < 1 (with flp < 0). However, we want the Skyrme functional to 
account for the finite-ranged pion, so it should not be equivalent to a short- 
distance expansion. Thus, the close correspondence suggests that the Skyrme 
functional is lacking and should be generalized. 

It is useful to compare results from the p only functional compared to the 
p and T functional, as in Fig. [TSland these tables: 

Op = Os E/A \/{r^) ap = 2as E/A ^/(r^) 
p Ym 2:87 ~p 03 3l0 

pT 7.65 2.87 pr 8.30 3.09 

We see very little difference in the Kohn-Sham observables, which are the bind- 
ing energy and the density distribution. However, the single-particle Kohn- 
Sham spectrum, which is not an observable, shows significant differences, as 
evidenced by Fig.[T9l (Note: The effective mass M* is closely related to single- 
particle levels.) We can show for a uniform system that the HP single-particle 
levels satisfy 

~ € - ^[(^ - ^ylrs + 2[u + l)al] , (105) 

and the pr result is the one corresponding to the spectrum from the full 
Hartree-Fock propagator. So the issue becomes how the full G (as shown below 
with the self-energy) is related to the Kohn-Sham Gks [M]; the closer they 
are, the better approximation Gks will provide for single-particle properties. 
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ii + 




To explore this connection, we add a non-local source £,{x',x) coupled to 
ip{x)il;^ {x') [we're back in Minkowski space here for no particular reason!]: 

Z[J ^] = e'^["''*l = j DijjDij}^ e' ^'^^'^ [c + J(x)^'' (x)^(x)+ j d^x' i,(x)i(x,x')i,'i (x')] 

(106) 

Writing r[p, i] = To [p, i] + Tint [p, , 



G{x, x ) — 



SW 



- ^ 



Gi,s{x,x') + G 



which is represented diagrammatically as: 



ks 



1 SFjnt ^ SFjnt 

i (5Gks 



Sp 



Gks , (107) 



I 

G 



Gv 



+ 



i.G 



ks 



Gv 



Gv 



Fig. 20. Full Green's function G in terms of the Kohn-Sham Green's function Gks. 



Now G and Gks yield the same density by construction; that is, yOks(x) = 
— zz^G^g(x, x'^ ) equals p(x) = — wG(x,x+). Here is a simple diagrammatic 
demonstration (the double line is minus the inverse of a single ph ring): 

But other single-particle properties (e.g., the spectrum) are generally different, 
since the last two terms in Fig. [20] will not cancel. 

We can ask whether the Kohn-Sham basis is a useful one for G. Or, more 
simply, ask how close is Gks to G. We find that it depends on what sources 
are used, as shown by the comparison of single-particle spectra in (jlOSp , with 
more sources implying less difference. This is a topic that merits further in- 
vestigation. 

The comparison of Kohn-Sham DFT and "mean-field" models often leads 
to misunderstandings, as when considering "occupation numbers" , because of 
a confusion between G and Gks- Figure W\\ suggests that occupation numbers 
are equal to or 1 if and only if correlations are not included. The Kohn-Sham 
propagator always has a "mean-field" structure, which means that (in the 
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Fig. 21. Schematic momentum occupation number n{k) for mean-field (Hartree- 
Fock) and with correlations. 



absence of pairing) the Kohn-Sham occupation numbers are always or 1. But 
correlations are certainly included in r[p\l (In principle, all correlations can be 
included; in practice, certain types like long-range particle-hole correlations 
may be largely omitted.) Further, n(k) = {<4i'^u) resolution dependent 
(not an observable!); the operator related to experiment is more complicated. 
Additional discussion on these issues can be found in [58] , 



3.1 Pairing in Kohn-Sham DFT 

There is abundant evidence for pairing in nuclei. The semi-empirical mass 
formula reproduces nuclear masses only with a pairing term (the last one): 



B{N,Z) = (15.6 MeV) 



1 - 1.5 



A 



A - {n.2Me\)A^/^ 



- (0.70 MeV)^ + (6MeV)[(-l)^ + {-l)^]/A^'^ , (108) 

which implies an odd-even staggering of binding energies (left panel of Fig.[22|). 
Other evidence is the energy gap in the spectra of deformed nuclei, low-lying 
2+ states in even nuclei (right panel of Fig. W^ . and deformations and mo- 
ments of inertia (the theory requires pairing to reproduce data) . 

Pairing is an example of spontaneous symmetry breaking, which is natu- 
rally accommodated in an effective action framework. For example, consider 
testing for zero-field magnetization M in a spin system by introducing an 
external field H to break the rotational symmetry. Legendre transform the 
Helmholtz free energy F{H): 

invert M ^ -dF{H)/dH =^ r[M] =^ F[H{M)] + MH{M) . (109) 

Since H = dF/dM — > 0, we look for the stationary points of F to identify 
possible ground states, including whether the symmetry broken state is lowest. 

For pairing, the broken symmetry is a C/(l) [phase] symmetry. The text- 
book effective action treatment in condensed matter is to introduce a contact 
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interaction [421143] : g tp'' tp'^ ipip , and perform a Hubbard-Stratonovich transfor- 
mation with an auxiliary pairing field A{x) coupled to -0^-0^, which eliminates 
the contact interaction. Then one constructs the IPI effective action r[A\ with 
A = (A), and looks for values for which Sr/6A = 0. To leading order in the 
loop expansion (mean field), this yields the BCS weak-coupling gap equation 
with gap A. 

The natural alternative here is to combine an expansion (e.g., EFT power 
counting) and the inversion method for effective actions [44j |45j |46j . Thus 
we introduce another external current j{x), which is coupled to the fermion 
pair density in order to explicitly breaks the phase symmetry. This is a nat- 
ural generalization of Kohn-Sham DFT [59l [60l [6T] . cf. DFT with nonlocal 
source [5^ 

So we consider a local composite effective action with pairing [64]. The 
generating functional has sources J,j coupled to the corresponding densities: 

Z[J,j] = £-^["''^1 = J D{%P^%P)e~^'^^'' [c+J{x)■4,li>^+3{x){^■^l+■4,^4,^)] _ ^-^-^^^ 

Densities are found by functional derivatives with respect to J and j: 

p(x)^(Vt(x)0(x))^., = , (111) 



and 



(j){x) = (4)\{x)^l{x) + i^i{x)^^{x)}jj 
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SW[J,j] 



6j{x) 



(112) 



The effective action r[p, 4>] follows as before by functional Legendre transfor- 
mation: 

r[p, = W[J,j] - J d^x J{x)p{x) - J d^xj{x)cj>{x) , (113) 

and is proportional to the (free) energy functional E[p,<p]] at finite temper- 
ature, the proportionality constant is /3. The sources are given by functional 
derivatives wrt p and <j): 

— -^=J(x and _iL^=j(x). 114 
(5p(x) d0(x) 

But the sources arc zero in the ground state, so we determine the ground-state 
p(x) and (/>(x) by stationarity: 



5E[p, ^] 



= 0. (115) 



This is Hohenberg-Kohn DFT extended to pairing! 

We need a method to carry out the Legendre transforms to get Kohn-Sham 

DFT; an obvious choice is to apply the Kohn-Sham inversion method again, 
with order-by-order matching in the counting parameter A. Once again, 

diagrams =^ W[JJ, A] = Wo[JJ] + XW,[JJ] + X^W^iJJ] + ■■■ 

assume =^ J[p, 4>, A] = Jo[p, + XJi[p, (j)] + X^J2[p, H 

assume =4> j [p, (j), A] = jo [p, 4>] + ^Ji [p, <P] + X^j2 [p, H 

derive =^ r[p, cj), A] = /q [p, </.] + AA [p, 0] + A^Fs [p, (/)] + ■■■ 

Start with the exact expressions for F and p 

r[p, = w[j,j] -jjp- jjcp, (116) 

and plug in the expansions, with p, <j) treated as order unity. Zeroth order is 
the Kohn-Sham system with potentials Jo(x) and jo(x), 

r^ip, 0] = Wo[Jo,io] - jjoP- jk<t>, (118) 

so the exact densities p(x) and ^(x) are by construction 
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(5Jo(a;) ' <5jo(a;) 
Now introduce single-particle orbitals and solve 



^o(x)-Mo io(x) 
jo(x) -/io(x)+^o 



Mi(x) 
fi(x) 



Wi(x) 
Wi(x) 



where 



V2 

^O(X) = - — + T4rap(x) - Jo(x) 



(119) 



(120) 



(121) 



This is just like Skyrme Hartree-Fock Bogliubov approach [8]. 

The diagrammatic expansion of the WiS is the same as without pairing, 
except now lines in diagrams are KS Nambu-Gor'kov Green's functions, 



G = 



(r,^|(x)V'j(x'))o {Tri^lix)4'i{x'))o i " 



F, 



ot 



(122) 



The extra diagram shown follows from the inversion (here it removes anoma- 
lous diagrams). In frequency space, the Kohn-Sham Green's functions are 



GL(x,x';c^)=^ 



i^,",(x,x';^) = -^ 



Wj(x)m*(x') t;j(x')w*(x) 



iu} — Ej iuj + Ej 



iuj — Ei 



ioj + Ej 



(123) 



(124) 



The Kohn-Sham self-consistency procedure involves the same iterations as 
in phenomenological Skyrme HF (or relativistic mean-field) when pairing is 
included. In terms of the orbitals, the fermion density is 



p(x) = 2^K(x)p, 



(125) 



and the pair density is (warning: this is unrenormalized!) 

^i^) = + ^»(x)<(x)] . 



(126) 



The chemical potential fig is fixed by J p{x) = A. Diagrams for r[p, <j)] oc 
Eq[p, (p] + £'int[p, (j)] yield the Kohn-sham potentials 



Jo(x) 



and jo(x) 



SEint [P,4>] 



(127) 
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3.2 Renormalization of Pairing 

When we carry out the DFT pairing calculation for a uniform dilute Fermi 
system, we find divergences almost immediately. The generating functional 
with constant sources /i and j is: 



^-wl^L,j] ^ 



J D{4,H) expiry d^x [^t (|. _ X! _ (128) 
■ ^°.^\^\^^^^+ (129) 



2 



(cf. adding an integration over an auxiliary field jD{A*\A)e i"^oi J ' , then 
shifting variables to eliminate for Z\*?/'|?/'|). There are new diver- 

gences because of j, e.g., expand W to 0{p): 



which has the same linear divergence as in 2-to-2 scattering. To renormalize, 
we add the counterterm ^CliP to C (see [33]), which is additive to W (cf. 
|Z\p), so there is no effect on scattering. 

We'll use dimensional regularization again, but generalize from DR/MS (as 
used by Papenbrock and Bertsch [65]) to DR/PDS, which generates explicit 
A dependence to "check" renormalization (by verifying that A dependence 
cancels). The basic free-space integral in D spatial dimensions is 



- — (A + ip), (130) 



where / 4!- ^ Renormalizing and matching free-space scattering yields 
for Co{A): 

Note: we recover DR/MS by taking yl = 0. As an exercise, you can verify that 
NLO renormalization in free space (left): 





implies that the corresponding sum of diagrams at finite density (right) is 
independent of A. 
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Now consider the Kohn-Sham noninteracting system for a uniform system, 
where we have constant chemical potential fiQ and pairing source jo (rather 
than spatially dependent sources). The bare density p is: 

1 dWo[fio,jo] _ 2 ■^^,2_ f 4-f^o \ (.oo\ 

and the bare pair density is; 

1 9M^o[Mo,jo] 2 ^ [ (fk jn 



PV djo ^ ^ J (2^)3 Ek 



In these expressions, jo plays role of constant gap; e.g., the spectrum is 

2M 

(See also Fig. [Ml) The divergence in (ps is illustrated in Fig. 

The basic DR/PDS integral in D dimensions, with x = ja/^o, is 



E, = ^(el~^or+fo , 4-^- (134) 



MA u( 
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We can check that in the KS density equation the A dependence cancels: 



1 dWol 



(fik 



1-4^ + ^ 



The KS equation for the pair density (j) fixes C^'^-' 
^ 1 dWo[] f d^k Jo 



0-/(l)+/io/(0) . (136) 



f3V djo 



(2^)3 Ek 



+ C'"'jo^-Jo/(0) + C(°)jo (137) 



so that 



MA 



(138) 



Calculating to rfi^ order, we find A<i<n[/0, 0] by first constructing all of 
the W"i<i<n[Aio(p, 0), jo(p, 0)], including additional Feynnian rules [46], 



So the procedure is to calculate /^i, ji from Fi, then use X^iLo-?* i ^ 
to find Jo- The renormalization conditions mean that there is no freedom 
in choosing Cq(A), so the yl's must cancel! In leading order, diagrams for 
A [p, (j)] = Wi [/^o(p, (^i), jo (/o, 4>)\ are 



Ej.tij.i'j, Ea.' lift',/ (i^j jo'i^'/j 



X— -{g)— --X 



and we choose and C^"-* to convert 0^ to the renormalized (j), yielding 



(139) 



The Fi dependence on p and 4> is explicit, so it is easy to find /ii and ji. 



1 ^^^ 1 (1) 



and 



1 dFi 



1 



9p 2 " (3V d(f) 

The "gap" equation then follows from j = Jq + ji = 0: 

d^k 1 



(140) 



(2-)^ v/(^r^7M^ 



(141) 

DR/PDS reproduces the Papenbrock/Bertsch result |65j (with x = Ijo/MoI) 

-1 



1 



-V2M^|as|(l + 2:2)'/'Pi" 



kpag 



/2 

4-61og2 2, 
1 — log a; 

TT TT 



(142) 
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Fig. 26. Contributions to the NLO energy density. 



and if kpUs < 1, then jo/a^o = (8/e^)e^^/^'''^l°='l holds to very good approxi- 
mation. 

The renormaUzed effective action F = Fq + Fi is 

^F = Ji4-^,,^E,) + lc^"^f,+^.op-Jo^+\cf^p' + \do'U'■ (143) 

We check for yl's again, 

^r = 0- I{2) + 2mo/(1) - {pi + jDm + + ■■■ (144) 

and find they do cancel: 
MA 



27r 



(-Mo(l - Jo/2/ii^) + - Mo - .70 + ^Jo) = . (145) 



To find the energy density, evaluate F at the stationary point jo = — ^|Cq"^^ \4> 
with /iQ fixed by the equation for p, yielding the same results as Papen- 
brock/Bertsch (plus an HF term) [64] . 

Life gets more complicated at Next-to-Leading Order (NLO), where de- 
pendence of I2 on p and (p is no longer explicit and analytic formulas for DR 
integrals not available. I2 at NLO is [Blj 



ici^ /tS 



and 



(27r)3 J (2^)3 J (2^)3 Ep + Ek+ Ep-q + Ek+q 

+ {uv)p {uv)k {uv)p^q {uv)k+q ] (146) 

-i^oyJ^^JpiukVk f + \^b(uI - vl)\' . (147) 
The UV divergences can be identified from 

2 _ 1 ik \ k^oo j^P 2 _ 1 /^i ^ a \ fc-oo ^ jIM^ 



(148) 
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and 



UkVk 



jo fc-»oo jpM 

'2Ek ^ F 



(2) _ 



For the renormalization at NLO to work, the bowtie diagram with the Cg 
{4:Tra^/M)A vertex must precisely cancel the A dependence from the beachball 
with Cq^'' = Anas/M vertices: 




+ 




A's cancel 



E UkVk 




1) 



(Note that the (^zj^"* vertex takes (j)B - 
and evaluate renormalized results without analytic formulas? 

Before addressing that issue, we first see how the standard induced in- 
teraction result [nni EZ] is recovered here. As jo 0, utVk peaks at /xq (see 
Fig. [211). At leading order (for T = 0), 



A's cancel 

E Uk'Vk' 

) How do we see cancellation of yl's 



^l/Af(0)|Co| 



(150) 



where we make the association jo ^ ^ lo ■ At NLO the exponent is modified, 
which changes the prefactor, A^lo ^ ^Lo/(4e)^/'^, using 



ri+r2 



+ 




1 



ii+i2 = ^\Co\ [1 - |Col(no)|k|=|k'|=fcFj ^ 



Further details can be found in [64j . An unexplored question is how does the 
Kohn-Sham gap compare to "real" gap? 

Now we return to the question of renormalizing in practice; an alter- 
native approach is to use subtractions. The NLO integrals with Ek — 
\J {tk — /io)^ + jo f^^c intractable, but we directly obtain a renormalized result 
with the substitution 



1 



E\ + E2 + £"3 + E^ 



1 



V 



Ei+E2+E^ + E4 e'l + €"2-4-4 
plus a DR/PDS integral that is proportional to A. When applied at LO 

MA 





■ 1 


V- 


firl 







2tt 



(151) 

I 

(152) 
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This is the same sort of subtraction used to ehminate Cq in the gap equation, 



M 
47ra„ 



1 

1^ 



(Pk 1 



(27r)3 el 



M 
47ra., 





■ 1 1 " 


/ (27r)3 





Any equivalent subtraction also works, e.g.. 



d^k 



d^k 1 



(2^)^ 4-^^o " 

So how do we renormalize the divergent pair (anomalous) density, 
= X! [< W^^i(x) + Wj(x)u*(x)] — > oo , 



(153) 



(154) 



(155) 



in a finite system? (Cf. the scalar density ps = ^(x)V-'(x) for relativistic 
Hartree theory). Answer: use the subtracted expression for cf) in the uniform 
system, 



'^^ d^k 



1 



1 



■Jo 



kc — >-oo 



(2-)^ VV(^".-Mo)^+Jo 



finite. 



(156) 



and apply this in a local density approximation (Thomas- Fermi): 



(x) = 2^ M,(x)v,(x) - jo(x)^^i^|^ with 



fcg(x) 
2A/ 



+ J(x) - ^0 



(157) 

This procedure was worked out by Bulgac and collaborators [60j EH [6T] . Con- 
vergence is very slow as the energy cutoff is increased, so Bulgac/Yu devised 
a different subtraction. 



1 



■Jo 



kc — *-oo 



finite 



(158) 



A comparison of convergence in uniform system for the two subtraction 
schemes ([T55)) and ([T55)) : 



10 

10-' 
^8 10-- 
-f 10- 

^ 10-' 
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i- ^lo=l■ 

Jo = 0-1 
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k 

subtract PHe^-iig) 
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shows dramatic iniprovement for the Bulgac/Yu subtraction. Bulgac et al. 
have demonstrated that this works in finite systems, and there have been 
recent practical apphcations. 

We finish this lecture with a brief mention of alternatives to a local Kohn- 
Sham formalism for pairing. One alternative is to couple a source to the non- 
local pair field [62l ES] : 

H — > H- Jdxdx' [D*{x,x')ip^{x)i:i{x') +Ii.c.] , (159) 

which yields essentially a two-particle-irreducible (2PI) effective action r[p, A] 
with A{x,x') — {ipi{x)ipi(x')) . Or one could use auxiliary fields: introduce 
A*{x)ip{x)ip{x) + H.c. via a Hubbard-Stratonovich transformation to obtain 
a IPI effective action in A{x) — {A{x)). By adopting a special saddle point 
evaluation, one can obtain Kohn-Sham DFT. Finally there is the possibility 
of deriving a density functional (without Kohn-Sham orbits) by direct renor- 
malization group evolution 41j. 

4 Loose Ends and Challenges plus Cold Atoms 

In this final lecture, we touch on some loose ends raised in previous lectures, 
and outline some of the plans and challenges for moving forward toward a 
microscopic DFT for nuclei based on effective field theory and renormalization 
group ideas and methods [55j. We'll also briefly consider cold atom physics 
and some recent work on density functionals for that problem. 

4.1 Tovirard a Microscopic Nuclear DFT 

We have outlined a framework for generating density functional theory based 
on effective actions. A key ingredient is a tractable hierarchy of many-body 
approximations to which we can apply the inversion method. A scenario for 
carrying this out has emerged, which combines chiral effective field theory 
(EFT) interactions with renormalization group techniques. While many chal- 
lenges remain, it is a plausible and systematically improvable path to a mi- 
croscopic nuclear DFT. 

This scenario goes like this: 

1. Construct a chiral EFT to a given order, including all many-body forces. 
At present, the NN chiral EFT has been worked out to N^LO [69l [70]. 
while three-body forces at the N^LO level are used. The latter will soon be 
extended to N'^LO and already the leading four-body force (which appears 
at N^LO) has been tested. To minimize the truncation error following 
Lepage's prescription, one should increase the cutoff regulator A until the 
truncation error is minimized. (Note: it is still a matter of investigation 
where the breakdown scale actually lies.) 
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2. Evolve the Hamiltonian to lower A with renormalization group (RG) 
methods. There are choices here, including the T^owfc approach, the Simi- 
larity Renormalization Group (SRG) [71] (75], and possibly simply a direct 
construction of the chiral EFT at a lower cutoff [73] . Cutoffs in the range 
of yl ~ 2 fm^^ appear to be appropriate for ordinary nuclei. One needs the 
consistent evolution of all interactions and other operators. As discussed 
in the first lecture, by decoupling high and low momentum the nuclear 
many-body problem becomes perturbative in the particle-particle channel, 
in stark contrast to the situation with conventional interactions [16J [75] . 

3. Generate the density functional in effective action form. A by-product 
of evolving to low momentum is that the convergence of the many-body 
diagrams no longer is critically dependent on the choice of single-particle 
potential. This opens the door to choosing it to maintain the density as 
in a Kohn-Sham approach. In the short term, a direct construction of the 
functional in the Skyrme form is possible via an adaption (and extension) 
of Negele and Vautherin's density matrix expansion (DME) [74l[75]. In 
the long term, chain-rule constructions will allow non-local effects to be 
included [see after ([88)1 ]. 

This program is well underway and is part of a larger project to construct and 
constrain a universal nuclear energy density functional (UNEDF). A detailed 
overview and an explanation of the DME will be available in a forthcoming 
publication [76] . 

We'll briefly describe the idea of the DME [74l[75], which starts by ex- 
pressing the Hartree-Fock energy using the density matrix. Recall that we 
take the best single Slater determinant in a variational sense 

Itf^Hp) =det{^,(x),i = x=(r,a,T), (160) 

to find the Hartree-Fock energy (suppressing a, r) : 




We can trivially express this in terms of the single-particle density matrix: 

p(ri,r2) = ^ 7/'l(ri)i/;a(r2) . (163) 

The idea is to write this in the Kohn-Sham basis (i.e., the ipa's are Kohn-Sham 
orbitals), so that it is compatible with the DFT diagrammatic expansion. If 
we change to R = i(ri -I- r2) and s = ri — r2, we can expand in s 



p(R + s/2, R - s/2) p(ri, ra)!,^, 



(164) 
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Fig. 27. Long-range effects contributing to energy density functionals. 

Negele and Vautherin obtained an expansion in terms of the fermion, kinetic 
energy, and other densities: 

r.) . «^p(R) + (i VV(R) - .(R) + |<cV(R) + ■ ^ ^) , 

(165) 

which leads to functionals of these densities, for which we can take the 
(5/5jo(R), 5/5t(R), etc. derivatives directly. (See also DME appHed to ChPT in 
nuclear medium by Kaiser and collaborators [771 [751 EH]-) This is clear at the 
Hartree-Fock level, but generalizations are needed for higher-order diagrams. 
These are also in progress [80] . 

There are some important open questions for this approach or any DFT 
treatment of finite nuclei. These include: 

• For pairing, the energy interpretation, number projection, renormalization 
in finite systems, and efficient numerical implementation. Also, a unified 
microscopic treatment of particle-particle and particle- hole physics. 

• DFT for self-bound systems. Self-bound systems have no external poten- 
tial, which implies that the true ground-state density is uniform! More gen- 
erally, how do we deal with symmetry breaking (translational, rotational 
invariance, particle number) and restoration. There has been little or no 
guidance from Coulomb DFT. There are analogous issues and methods 
for effective actions, namely soliton zero modes and projection methods. 
Work on an energy functional for the intrinsic density is in its infancy ^81J . 

• Long-range effects, as illustrated schematically in Fig. [27l This includes 
long-range forces (i.e., pion and Coulomb exchange) but also long-range 
correlations. The latter can be understood as non-localities from near-on- 
shell particle-hole excitations, as in the lower diagrams pictured in Fig. 1271 
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4.2 Covariant DFT 

Thus far, our discussion has included only nonrelativistic EFT and DFT. 
However, there is a successful phenomenology based on relativistic mean field 
energy functional [121 [83l [84] . Can we make a connection? 

In principle we could proceed by deriving a covariant EFT. We start by 
observing that all low-energy effective theories have incorrect UV behavior. 
Sensitivity to short-distance physics is signalled by divergences but finiteness 
(e.g., with cutoff) doesn't mean there is not sensitivity! One must absorb 
(and correct) sensitivity by renormalization. Instances of UV divergences for 
low-energy nuclear physics are 

nonrelativistic covariant 
scattering scattering 
pairing pairing 

anti-nucleons 




Thus, there is an additional source of divergences in the covariant case from 
the "Dirac sea" . 

Gasser, Sainio, Svarc [85] derived chiral perturbation theory (ChPT) for 
ttN physics using relativistic nucleon degrees of freedom. But they found that 
the loop and momentum expansions no longer agree (as they do in nonrela- 
tivistic ChPT), which means that systematic power counting was lost. The 
heavy-baryon EFT restores power counting by a 1/M expansion, and has been 
the basis for nonrelativistic NN EFT treatments "86]. 

However, Hua-Bin Tang [Ff] (and with Paul Ellis [881) observed: 

". . . eft's permit useful low-energy expansions only if we absorb all 
of the hard-momentum effects into the parameters of the Lagrangian." 
"When we include the nucleons relativistically, the anti-nucleon con- 
tributions are also hard-momentum effects." 

They advocated moving the "Dirac sea" physics into the coefficients, thereby 
absorbing the "hard" part of a diagram into parameters, while the remaining 
"soft" part satisfies chiral power counting. The original irN prescription by 
Tang and Ellis (expand, integrate term-by-term, and resum propagators) was 
systematized for ttN by Becher and Leutwyler under the name "infrared reg- 
ularization" or IR [89]. It is not unique; e.g., Fuchs et al. have used additional 
finite subtractions in DR [SU] . The extension of IR to multiple heavy particles 
by Lehmann and Prezeau [91j . with a convenient reformulation by Schindler, 
Gegelia, and Scherer [9^, offers the possibility of a working covariant EFT. 
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If we restrict our attention to purely short-ranged, natural interactions, 
there are tremendous simplifications. In particular, tadpoles and NN loops 
in free space vanish! For example, leading order (LO) has scalar, vector, etc. 
vertices, 

which we designate as 



and consider all possible diagrams at NLO: 




Only the particle-particle loop diagram survives IR and all of the others pic- 
tured here vanish. Since only forward-going nucleons contribute in the end, 
one obtains the same scattering amplitude as in nonrelativistic DR/MS for 
small k. 

Unlike QED DFT, "no sea" for nuclear structure is a misnomer; one should 
include NN "vacuum physics" in coefficients via renormalization. But note 
that requiring renormalizability at the hadronic level corresponds to making 
a model for the short-distance behavior, which has proven to be a poor model 
phenomenologically. Fixing short-distance behavior is not the same thing as 
throwing away negative-energy states. For a long time, people searched for 
unique "relativistic effects" ; these were largely misguided efforts. 

The further investigation of covariant EFT and its extension to DFT is 
motivated by the successes of relativistic mean-field phenomenology and other 
arguments about low-energy QCD. But there is much to be done for it to be 
competitive with the nonrelativistic EFT. 

4.3 DFT for Cold Atoms with Large Scattering Length 

Finally we return to the large scattering length problem, which is realizable 
with cold atoms. The total cross section for scattering is expressed in term of 
partial-wave phase shifts as 



1=0 



(167) 
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Recall that an attractive potential pulls the asymptotic wave function (outside 

the potential) in, by an amount at each energy or k called the phase shift. At 
low energy (A = 27r/fc 1/-^)) the S-wave phase shift ^o(^) satisfies: 



where ao is the "scattering length" and Tq is the "effective range" . The effective 
range expansion for low-energy scattering goes back to Schwinger and Bethe 
and others. The effective range expansion typifies the general principles we 
have stated for EFT's: If a complicated potential produces scattering with 
a given ao and ro, we can replace it by a simpler potential with the same 
values and everything agrees at low energies. In general, the effective-range 
expansion is reproduced and extended by EFT. 

Having a bound-state or near-bound state at zero energy means large 
scattering lengths (ao — > ±00). For kR — > 0, the total cross section is 



We are particularly interested in cases where there is a bound-state near zero 
energy or there just misses being a bound state. These pictures: 



are a reminder of the interpretation of the scattering length in terms of the 
intercept of the zero energy scattering wave function, which is a straight line 
outside the potential. For potentials that just have a bound state, the wave 
function just turns over and ao is large and positive. If the potential just fails 
to have a bound state, it doesn't quite make it to horizontal and ao is large 
and negative. 

At low energies, depending on the size of k times ao, the cross section first 
goes like the square of ao and then saturates at the unitarity limit. So if we 
could adjust the depth of the bound state, we can control the "strength" of 
the interaction in a sense. This is possible for atoms by changing an external 
magnetic field to produce resonant scattering. For QCD is it possible to adjust 
the quark mass theoretically so that m„ changes and the nuclear ao can be 
tuned to ±00! 



kcot6o{k) — — + Irrok"^ + ■ . ■ 



(168) 




(169) 
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Let's consider the large scattering length many-body problem, which 
means we have an attractive two-body potential with oq — s- oo. If i? <C l/fcp ^ 
Oq, as in this figure, 




then we expect scale invariance (since we lose both R and ao as possible scales). 
This means that the energy and superfluid gap should be pure numbers times 

TP 3 

-'^FG - 52M- 

Recall that for the natural scattering length case, EFT power counting led 
to an organized perturbative expansion: 



P/2 + k 



P/2 + k' 




P/2-k P/2-k' 

iT{k, cosB) 






with Co = ^'^0 and fcap ^ 1. But in the large scattering length limit, 
fcao ^ 1 so the bubble series diverges. This is not a difficulty in free space, 
because the geometric sum of bubbles is easily performed. This sum yields the 
fo{k) expansion by keeping fcao to all orders and expanding the rest: 





foik) cx 



-1 



— l/ao + rok'^/2 — ik l/a^ + ik 



1 



ro/2 



l/ttQ + ik 



(170) 



With a natural oq and a perturbative expansion, we found the DR/MS 
(minimal subtraction) scheme particularly convenient. With large ao, we need 
a new renormalization scheme. DR/PDS was proposed by Kaplan, Savage, 
Wise [93| and counts /Lt ~ fc: 
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Ea 



Pairing gap (A) = 0.9 E 






Fig. 28. GFMC results from Chang et al. [94] for the unitary fermion system. 




Co(p) 



(171) 




(27r)3 P 



An 



In medium, each additional Cq vertex gives a factor 



CoikF) 



(172) 



(173) 



which means that all Co diagrams are leading order! Thus, we are told to sum 
all many-body diagrams with Cq vertices. This is only possible numerically 
(or possibly with an additional expansion). 

So we turn to numerical calculations. GFMC results from Chang et al. 
in which one extrapolates to large numbers of fermions, are shown in Fig. 1281 
They find the energy per particle is E/N = 0.44(l)i?FG- Diffusion Monte 
Carlo (DMC) results [95] , with a square- well potential tuned to ao ^ oo and 
an extrapolation to large numbers of fermions, find a similar result, namely 
an energy per particle of E/N = 0.42(l)i?FG- 

Thomas Papenbrock has considered DFT for the unitary regime [9^ . He 
assumes a simple, constrained form of the density functional. 



m - - 

m 



N 



2m 



off 



cp 



5/3 



+ ^muj^r^p , (174) 



with non- localities and gradient terms via the effective mass rrics- The pa- 
rameters in E[p] = J dx£[p{x.)] can be fit for = 2 to exact results for two 
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Fig. 29. Papenbrock results from a density functional for the unitary regime. 



fermions in harmonic trap [Busch et al (1998)] 

1 I ^^^lAl'i 



I = ^/h/muj , E — 2Tiu! , 



and a gaussian center-of-mass wave function. The result is 



Pcxact(^) 



4 i^-2(r/lf 



^3/2/3 



dx ( 



(175) 



(176) 



Results of such fits are shown in Fig. [23 They predict <^ = 0.42 and 
m/mes = 0.69 from the best fit. The value for ^ is amazingly close to the 
Monte Carlo result. Papenbrock and Bhattacharyya [HTKSH] consider correc- 
tions for an LDA density functional close to the unitary limit 



Cl 



ap 



1/3 



1/3 



(177) 



Again, this is fit to the harmonically trapped two-fermion system. Results are 
given in Fig.[30]and show impressive agreement with Monte Carlo calculations. 

Finally, there are interesting investigations of the constraints of general 
coordinate and conformal invariance by Dam Son and collaborators |99i ilOO] • 
They ask: Is there more than scale invariance for the unitary Fermi gas? The 
symmetries can be exposed by adding a background gauge field and curved 
space with metric (7y (t,x): 



S 



dt dx 



-V-ta^V - |-(5. + iA,)^^{d, - lA,)^ (178) 
2 2m 



2 



2r2 



(179) 
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Fig. 30. Papenbrock/Bhattacharyya DFT results for finite scattering lengtii a com- 
pared to Monte Carlo calculations. 



This is more than scale and Galilean invariance! Direct consequences include 
extra constraints on Ccs at NLO, which naively involve five arbitrary func- 
tions but these symmetries show there are only three. For the unitary Fermi 
gas, three constants from scale invariance are reduced to two constants from 
conformal invariance. This leads us to ask: What additional constraints can 
we find for the energy functional? This and the other open questions are ripe 
for investigation! 
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